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ABSTRACT 

This Letter investigates the origin of the asymmetric magnetic field line geometry in 
the ergospheric disk (and the corresponding asymmetric powerful jet) in 3-D perfect 
magnetohydrodynamic (MH D) numerical simula tions of a rapidly rotating black hole 
accretion system reported in 'Punsl v et all (|2QQ9f ). Understanding, why and how these 
unexpected asymmetric structures form is of practical interest because an ergospheric 
disk jet can boost the black hole driven jet power many-fold possibly resolving a 
fundamental disconn ect between t he energy flux estimates of po werful quasar jets and 
simulated jet power (|Punslvll2Qllf ). The new 3-D simulations of iBeckwith et aj (|2QQ9f ) 
that were r un with basically the same code that was used in the simulation discussed in 
iPunslv e~ al (2009) describe the "coronal mechanism" of accreting poliodal magnetic 
flux towards the event horizon. It was determined that reconnection in the inner 
accretion disk is a "necessary" component for this process. The coronal mechanism 
seems to naturally explain the asymmetric ergospheric disk field lines that were seen 
in the simulations. Using examples from the literature, it is discussed how apparently 
small changes in the reconnection geometry and rates can make enormous changes 
in the magnetospheric flux distribution and the resultant black hole driven jet power 
in a numerical simulation. Unfortunately, reconnection is a consequence of numerical 
diffusion and not a detailed (yet to be fully understood) physical mechanism in the 
existing suite of perfect MHD based numerical simulations. The implication is that 
there is presently great uncertainty in the flux distribution of astrophysical black hole 
magnetospheres and the resultant jet power. 
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The ultimate nature of the power output of a black hole 
magnetosphere is highly dependent on two poorly under- 
stood circumstances, the source of plasma injection on the 
magnetic field lines that thread the e vent horizon and the 
fate of accreted vertical magnetic flux (jPunslv and Coronitil 
Il990a.b). In the past eight years, the use of perfect mag- 
netohydrodynamic (MHD) numerical simulations have been 
developed in the scientific community to help understand 
these issues. In order to establish and maintain the event 
horizon magnetosphere, the simulations must rely on the 
numerical artifice of a mass floor, local mass injection to 
establish a minimum density. In spite of some discussion of 
anecdotal mass floor examples, perfect MHD simulations are 
not likely to resolve the first issue since a mass floor violates 
rest-mass and energy - momentum conservation and neces- 
sarily contradicts the perfect MHD assumption (jMcKinney, 



l2006bl ). However, recently MHD simulations have shed some 
light on the second point. The mystery of how magnetic flux 
was accreted to the b lack hole in MHD nu merical simula- 
tions was revealed in iBeckwith et all (|2009h through their 
meticulous high resolution numerical work. In this Letter, 
the insight provided by Beckwith et al (2009) is used to ex- 
plain the strange ergos pheric disk fie l d top ology seen in the 
simulations reported in IPunslv et all (|2009l) that resulted in 
powerful one-sided jets that dominate the total energy out- 
put in the jetted system with a location that changed hemi- 
spheres in different time snapshots (Punsly 2007a, b). 

The fate of accreted flux is perhaps the most critical is- 
sue in understanding the power source for black ho le driven 
jets. It was noted in IPunslv and Coronitil (|l990bl ). that if 
vertical, magnetic flux accretes, it is not clear where it ends 
up. It was argued that reconnection of vertical flux would be 
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determinant to the final magnetic field configuration since 
the black hole is effectively a sink with infinite capacity for 
mass, but with a very limited capacity to accept magnetic 
flux. One possible field configuration that could result from 
reconnection produced a disk in the ergosphere that can 
drive powerful jets. If this happens, the power source for the 
jet drastically increases in effic iency, so this configu r ation is 
of pr ofound interest in AGN (Nemmen et al' '20071: IPunslvl 
I2OIII ). However, there is much scientific uncertainty in the 
reconnection geometry and rate expected in an accretion 
flow near a rapidly rotating black hole, since our experimen- 
tal experience is based on very different environments, the 
solar corona, the solar wind, the Earth magnetosphere and 
magneto-tail, and magnetic confinement devices for ther- 
monuclear fusion. The situation was further complicated by 
the argument that the existence of coherent vertical flux 
within the dense accreting gas is inh ibited by the turbulent 
magnetic diffusivi ty of the plasma (|van BallegooijenI Il989l : 
iLubow et allllOQ^ l. Thus, it was not even clear if the notion 
of a large scale field associated with an accretion flow was 
viable. 



However, recent numerical simulations have shown that 
vertical flux accretion can occur, but not by diffusing 
through the disk, but by a two step process called the 
"coronal mechanism". The first step is the transport to- 
wards the black hole of an inwardly stretched, disk- anchored, 
poloidal loop through the low turbul ence, coronal la yer just 
above the disk, the "hairpin" field in lSeckwith et all ((20091 ) . 
T his coronal transport is similar to the mech anism proposed 
i n iRothstein and Lovelac^ (|2008l l. but see iBeckwith et all 
(|2009l ) for some dynamical differences. Then loops of mag- 
netic field in the inner accretion flow reconnect with the half 
of the "hairpin" at mid- latitudes, allowing it to contract into 
the black hole, leading to one sign of field threading the black 
hole (see section 2). Reconnection is apparently "necessary" 
for vertical flux accr etion in a ll the numer ical simulations 
of th is familv Beckwith et all (12009): Hawlev and KrolikI 
(|2006l l: iMcKinnev and BlandfordI (|2009l l. but the geometry 
of the rec onnection site is very differen t than what was de- 
scribed in IPunslv and Coronitil (Il990bl ). In spite of this geo- 
metric difference, an ergospheric disk does form in some sim- 
ulations. In the following, the coronal mechanism is shown to 
naturally explain the asymmetric magne tic field observed i n 
3-D simulations of the ergospheric disks (jPunslv et alll2009h . 



work derives from a family of simulations 
on a constrained transport MHD code on 
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the Kerr spacetime background that has been de- 
scribed numerous times in the literature, so this 
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is not reproduced here (De Villiers and Hawlev 2003a; 
liers et al 2003b, 2005: Hirose et al 2004: Krolik et al 
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20051 : iHawlev and KrolikI l2006l : iBeckwith et all l2008l). The 
particular simulation, KDJ, discus sed in | Punslv etiil 
(l2009h was described i n detail in iKrohk et all (12005 



iHawlev and Krohkl (|2006), which the reader should consult 
for particulars. KDJ simulates accretion onto a rapidly ro- 
tating black hole with an angular momentum per unit mass 
of a = 0.99 M in geometrized units. This Letter focuses on 
the very complicated twisted topologies that are involved in 
the reconnection events in the simulation. 



1 FIELD LINE TOPOLOGY IN THE 
ERGOSPHERIC DISK 

A surprising feature of 3-D MHD numerical simulations is 
that the ergospheric disk was threaded mainly by the type I 
field lines instead of type HI field lines in the nomenclature 
introduced in IPunslv et all (|2009l ) . In the left hand panel 
of Figure 1, the type I vertical magnetic field lines emerge 
from the inner equatorial accretion flow. The false color plot 
is a contour map of a 2-D cross-section of the density in 
Boyer-Lindquist coordinates expressed in code units. The in- 
terior of the inner calculation boundary (just outside of the 
event horizon) is grayish- white. These field lines are distin- 
guished by connecting to the Poynting jet in one hemisphere 
only, with the other end spiraling around within the accret- 
ing gas in the opposite hemisphere. Another distinguishing 
feature of a type I field line is that the azimuthal direc- 
tion of the magnetic field changes direction as the field line 
cr osses the midplane of the accretion flow. It was expected 
in Punsly and Coronitii (1990b) that the direct advection of 
magnetic flux would produce the type HI field topology de- 
picted in the right hand frame of Figu re 1. But the s e wer e 
rare in the 3-D simulation reported in IPunslv et all (|2009l ). 
Unlike the type I field lines, the type HI field lines connect 
to both sides of the bipolar jet. 

One should consult Figure 3 of IPunslv et all (|2009l l to 
appreciate the physical significance of the one sided type 
I field lines from the ergopsheric disk. These field lines di- 
rectly equate to a powerful jet of Poynting flux. When this 
jet f orms, it swamp s the power output from the event hori- 



zon (|Punslvll20073 lbh. In fact, it even suppress es the event 



horizon jet power to some degree 



3ven suppress e; 
(|Punslvll20nT ). 



2 RELEVANT ASPECTS OF THE CORONAL 
MECHANISM 

In this section, the coronal mechanism is reviewed and asym- 
metric variations in the process are noted that are relevant to 
the followi ng. The c oronal mechanism was detailed in Figure 
11 of Beckwith et al (2009). Som e important geometric sim- 
plifications were implemented in IBeckwith et all (|2009l l that 
greatly improved the clarity of presentation. In that paper, 
the primary focus was on detailed discussions of symmet- 
ric pairs of "hairpin" field lines, where one approaches from 
just above the accretion disk from the North and a matching 
"hairpin" field line approaches from the South. Furthermore, 
a non-rotating black hole was chosen so there is minimal az- 
imuthal twisting near the event horizon. Ostensibly, for the 
sake of simplicity of pr esentation, all t he 3-D data was av- 
eraged over azimuth in IBeckwith et all (|2009l l , so there were 
simple field line topologies like closed 2-D poloidal loops in 
the accretion disk (see their Figure 11). 

After the leading edge of the hairpins have penetrated 
the event horizon, reconnection occurs in the equatorial por- 
tion of the hairpin with a closed 2-D loop. This leads to an 
increase in the magnetic flux in the accretion flow vortex, 
or "funnel", in the polar region beyond the event horizon 
radius in both hemispheres. 

It is noted here that there are t wo intriguing aspects o f 
the field line evolution presented in IBeckwith etail (|2009l l 
that should be relevant to 3-D field line reconnection in the 
Kerr geometry. 
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Figure 1. Plots of ergospheric disk field lines in the simulation, KDJ, from Punsly et al (2009). For details of how the field lines were 
plotted and the false color 2-D density cross section, please refer to that paper. It was expected that the type III ergospheric disk field 
line topology (on the right) would dominate the vertical fiux through the ergospheric disk. Instead 3-D MHD numerical simulations found 
very few type III field lines, but many type I field lines (on the left). It should be noted that the images are 3-D. They were created 
in the 3-D visualization tool, Paraview 3.3, as 3-D images and the "camera" is placed in the equatorial plane (i.e., if the y-axis points 
toward the camera, the gas density is plotted in the x-z plane). The field lines that pass behind the x-z plane are made darker (i.e. the 
x-z plane is partially opaque, 75%) to make it easy to distinguish them from the magnetic coils in front of the x-z plane. However, near 
the poles they are not darker because a conical region of the numerical grid is excised in that region (therefore, there are no points in 
the x-z plane above the poles of the horizon and no opacity) as a numerical expedience as discussed in Hawley and Krolik (2006) and 
Punsly et al (2009). 



(i) There is not always symmetry in the hairpin accretion 
between the northern and southern hemispheres. A cluster 
of field lines in a "hairpin" configuration often approach 
from one hemisphere at a time. Th is is manifested in the on- 
line animation for Beckw ith et all (^2009), especially around 
23 seconds, where one sees an excess of organized flux in 
the funnel, in the southern hemisphere of the event hori- 
zon. This asyr ametry also appears t o a much lesser degree 
in Figure 11 of Beckwith et a] (|2009l ). Note that the anima- 
tion that accompanies the 2-D simulation around a rapidly 
spinning Kerr black hole in IBeckwith et all (|20Q8l ) also shows 
asymmetric hairpin accretion. This supports the notion that 
asymmetric hairpin field line accretion is not dependent on 
spin. This does not mean that reconnection will proceed sim- 
ilarly in 2-D and 3-D and independent of spin. Instead the 
simulations show only that the potential pre-reconnection 
field configuration of asymmetric accretion of hairpin field 
lines seems to occur in either 2-D or 3-D and for high spin 
and low spin. 

(ii) Some of the hairpins can actually be buoyant near the 
horizon and move away from the black hole, never penetrat- 
ing the horiz on (eg;., the hai r pin ^ 30° below the equator in 
Figure 11 of IBeckwith et all (|2009l ) moves outward between 
t= 14600 M and t= 14640 M) 



up into toroidal coils ("H irose et all2004l : IPunslv et a]|l2009l l. 

There are no simple poloidal loops in 3-D li ke there were in 
the az imuthally averaged data in Figure 11 of lBeckwith et aj 
{"2009). Since topology is critical to reconnection, the 2-D ex- 
pedience is not implemented here for the sake of accuracy 
and the expense of complexity. The analog of a 2-D poloidal 
loop in the Schwarzschild geometry in the 3-D Kerr (rotat- 
ing black hole) geometry would seem to be one of the many 
twisted coils that p ermeate the accretion flow iHirose et all 
(2004): IPunslv et al (2009). The "loop" from t = 9840 M (in 
geometrized units) of KDJ that is plotted in white in the left 
hand frame of Figure 2 is typical of most of the field lines in 
the inner accretion disk for high spin black holes. It spirals 
near the black hole then as it expands vertically, it leaves the 
disk and penetrates the corona. The twisted loop connects to 
large distances through the corona, presumably closing far 
away. There are variations of this topology in which one end 
of the coil stays in the disk while the other en d permeates 
the corona or event horizon (IPunslv et a]|l2009l ). 

There are no fine time resolution snapshots of the sim- 
ulation KDJ (the data is sampled every 80 M), so we can- 
not see the asymmetric type I field lines of Figure 1 form- 
ing. There is only circumstantial evidence as to the chain of 
events: 



3 RECONNECTION GEOMETRY IN THE 3-D 
ACCRETION NEAR RAPIDLY ROTATING 
BLACK HOLES 

The reconnection aspect of the "coronal mechanism" de- 
scribed in IBeckwith et all |2009l l relies on the topology of 
poloidal loops in the inner accretion disk. However, in the 
case of interest here, 3-D around a rapidly spinning black 
hole, almost all the flux in the inner accretion flow is twisted 



(i) Most of the disk field lines in KDJ are twisted coils 
near the event horizon that can have significant random in- 
flections due to turbulence and never l eave the disk or corona 
(|Punslv et a]|l2009l : lHirose et a]|l2004l ). 

(ii) In KDJ, large vertical flux bundles are entrapped in 
the inner regions of the accretio n flow tending to be far more 
pronounced in one hemisphere (I Punslv et alll2009l ). 

(iii) In KDJ, the ergospheric disk flux shows no tendency 
for a preferred hemisphere. They can be in either hemi- 
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Figure 2. A potential reconnection site in KDJ is displayed in the left hand frame. The red field line is the accreting, slightly twisted 
coronal hairpin. The white field line is a typical twisted structure that is the 3-D equivalent of a poloidal loop in a 2-D Schwarzschild 
representation. The inner calculational boundary is blue. The post reconnection topology is depicted in the right hand frame. The red 
curve is a buoyant hairpin. The white curve is a type I field line that threads the ergospheric disk as in Figure 1. 



sphere or bo th hemispheres in an individual time snapshot 
(jPunslv et a l 2009). 

(iv) The animations of simulations show that significant 
poloidal flux can accrete to the inner regions of the accretion 
flow thr ough the corona in the form of hairpins in an episodic 
fashion (Beckwit h erall2008l . lioool ). 

(v) These hairpins can arrive at the inner ed ge of the flow 
in as ymmetric north/south configurations (|Beckwith et all 
l2008l . 120091 

(vi) The topology of the poloidal flux near the event 
horizon was shown to be determined by reconnection with 
simple poloidal loop s in 2-D azimuthally averaged data 
(jBeckwith et alll2009l ). 

(vii) Figure 11 of iBeckwith et all (|2009l l not only shows 
the asymmetric accretion of hairpins to the event horizon, 
but also that some hairpins near the black hole can actually 
become buoyant and move away from the black hole. 

These seven facts can be used to consider the time evolution 
in the twisted geometry of KDJ that results in the type I 
field line topology in Figure 1. Figure 2 simulates the most 
plausible scenario based on these seven results above which 
is used as a surrogate for actual fine time scale data sam- 
pling. The left hand frame shows a typical twisted accretion 
disk coil in white being approached by a coronal hairpin 
in red near the black hole. Notice that the coronal hairpin 
is azimuthally twisted in KDJ. Reconnection is very com- 
plicated i n a twisted 3-D e nvironment and is not well un- 
derstood (|Pontin et alll20l"ll V However, the configuration as 
drawn forms a natural reconnection site (an X-point). We 
expect both types of field lines in Figure 2 to exist from 
points 1, 4 and 5 above. The elements required for the pre- 
reconnection geometry in the left hand frame of Figure 2, 
commonly occur in this family of simulations. Thus, it is 
reasonable to expect that these field configurations coexist 
in proximity at various times and these potential reconnec- 
tion sites should not be rare. However, the reconnection rate 
in such a complicated geometry that does not proceed by 
a physical mechanism, but through numerical diffusion, is 
very uncertain. By points 2, 3 and 6 above, reconnection 
must have occurred in KDJ as depicted in the right hand 
frame of Figure 2. The white curve in the right hand frame is 
poloidal flux through the equatorial plane of the ergosphere 
in one hemisphere in analogy to the left hand frame of Fig- 
ure 1 and the red curve would be a buoyant hairpin field line 
that moves out in the corona consistent with point 7. 



A significant difference between the topology resulting 
from the reconnection in Figure 2 compared to that in Fig- 
ure 11 of iBeckwith et all (|2009h , is that in Figure 2 the re- 
connection is happening before th e hairpin penetrates the 
event horizon and in Figure 11 of IBeckwith et all (|2009l l it 
occurs after the hairpin penetrates the horizon. This indi- 
cates that the coherent flux transport rate combined with 
the reconnection rate and twisted 3-D field line geometry 
(which aflPects the reconnection rate) might determine if a 
field line penetrates the event horizon or the inner accretion 
flow when and if reconnection occurs. The flnal fleld line 
topology depends on the balancing of reaction rates (recon- 
nection and transport) as well as internal dynamics (that 
affect fleld line shape) that are determined by the numerical 
simulation. 



4 DISCUSSION 

The Letter shows that the "coronal mechanism" for flux 
transport in simulations of black hole accretion provides a 
plausible explanation for the one sided ergospheric disk fleld 
lines in the high spin 3-D simulation KDJ. It therefore ex- 
plains the strange phenomenon observed in KDJ that the 
black hole driven jet Poynting flux was very one sided, jump- 
ing from side to side and emanating primarily from the er- 
gospheric disk. 

An otherwise almost identical simulation to KDJ that 
includes additional artiflcial diffusion terms in the equations 
of continuity, energy con servation, a n d mo mentum conser- 
vation (as described in iDe VilliersI (|2006l ) ) do not show 
these one sided ergospheric disk structures IPunslvl (|201l[ ) . 
This is seems to indicate a change in the the reconnec- 
tion process that is driven either directly or indirectly by 
the numerical diffusion. In support of this interpretation, 
the force-free simulations of an initially uniform fleld in 
iKomissarovl ()2004l ) show magnetic flux threading the ergo- 
spheric equatorial plane near the black hole, yet the same 
initial state that is time evolved in a different force-free 
code with a slower ansatz for the reconnection rate shows 
no magneti c flux threading the equatorial plane near the 
black hole (jMcKinnevI l2006al ) . The implication is that the 
global topology of the black hole magnetosphere is highly 
dependent on the time evolution driven by reconnection. 
This is not a trivial circumstance because in MHD simula- 
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tions reconnection is dependent on numerical diffusion. The 
situation is rendered even more ambiguous by the comph- 
cated twisted magnetic field line geometries in 3-D simula- 
tions around rapidly rotating black holes (eg. Figure 2). The 
complications of far less intricate 3-D field line t opology have 
been recog;nized in solar and plan etary physics (|Pontin et all 
I2OIII : IWilmot-Smith et all I20T0I ). Similarly, some detailed 
numerical modeling ha s shown the need for 3-D to p roperly 
descr ibe reconnection (|Kowal et all 1 20091 : iKulpa-Dvbel et all 
I2OIOI ). Furthermore, in high energy environment s radiation 
effects might also be crucial for a proper treatment lUzdenskvl 
(|201ll ). This might be relevant because radiation is not for- 
mally considered in these MHD codes. In summary, there 
are many potential sources of uncertainty in the reconnec- 
tion induced topology in these simulations. 

It should be noted that in principle, the choice of coor- 
dinates (eg. Boyer-Lindquist, as in the simulations discussed 
here, or Kerr-Schild) should be irrelevant to the results ob- 
tained. What is likely more relevant is the grid scale and 
sources of numerical diffusion in the solution methodology. 
A very basic consideration is that once a field line that is 
frozen-in to the plasma is accreted through the event horizon 
(technically, the inner calculational boundary that is just 
outside the event horizon in Boyer-Lindquist coordinates), 
it should be unable to migrate out of the horizon by moving 
outward into the equatorial plane if causality is maintained 
by the frozen-in plasma inside the horizon. Field lines can 
only be extracted from the horizon by local reconnection, 
not diffusion back into the equatorial plane. As such, a black 
hole saturated with magnetic flux is a nontrivial "boundary 
condition" for further magnetic flux accretion of the same 
orientation. For example what happens when there is re- 
peated injection of m agjnetic flux as in the simulations of 
llgum enshchevi (|2008h . The suite of simulations considered 
here have a finite amount of flux (either vertical or the lead- 
ing edge of accreting loops) that accretes towards the black 
hole and a boundary condition that flux cannot leave or 
enter through the outer boundary. However, if that bound- 
ary condition is changed to be one of continuous injection 
of flux then the magneti c topolo gy might change near the 
black hole. In particular. iBeckwith et af(j2009l ) suggest that 
the funnel and horizon field stre ngth is set by the gas and 
magnetic pressure in the disk and lTchekhovskov et ail (|2010l l 
claim that it is set by the ram pressure in the disk. Once 
this very finite value is achieved, what happens when more 
fiux is injected into the accretion stream and approaches 
the horizon? Does this preferentially move the reconnection 
sites out to near the mid plane of the ergosphere? This might 
raise the power le vel to be consistent with the more powerful 
radio loud AGN (jNemmen et alll2007l : |Punslvll201lh . 

In the absence of more robust treatments of recon- 
nection, this study suggests two new simulations that are 
achievable with the present numerical codes. One would be 
to rerun KDJ with the original code, but with high density 
time sampling for a portion of the run so that we can see how 
the ergospheric disk fiux gets established. It would also be 
an interesting simulation to explore the possibility of contin- 
uous vertical fiux injection from the outer boundary around 
rapidly rotating black holes to see if the pattern of reconnec- 
tion changes and if the magnetospheric distribution of fiux 
changes, possibly enhancing the ergospheric disk. In the fu- 
ture, more realistic reconnection in the presence of resistiv- 



ity ne eds to modeled near black holes as in IPalenzuela et al 
(|2009h . One might also find environments suitable for recon- 
nection in the ergospheres of other simulation geometries 
that can produce relativistic jets. For example, th e colli - 
sions of black holes as discussed in IPalenzuela et all (|2010h 
or the inner region s of a tilted accretion disk as modeled in 
iFragile et aj (|200/t ) are fertile areas to pursue. 
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